Kinetics of spin coherence of electrons in n-type InAs quantum wells under intense 

terahertz laser fields 



J. H. Jiang,i'2 M. W. Wu,i'2.E] and Y. Zhoui 

' Hefei National Laboratory for Physical Sciences at Microscale, 
University of Science and Technology of China, Hefei, Anhui, 230026, China 
Department of Physics, University of Science and Technology of China, Hefei, Anhui, 230026, China 

(Dated: August 20, 2008) 

Spin kinetics in n-type InAs quantum wells under intense terahertz laser fields is investigated by 
developing fully microscopic kinetic spin Bloch equations via the Floquet-Markov theory and the 
nonequilibrium Green's function approach, with all the relevant scattering, such as the electron- 
impurity, electron-phonon, and electron-electron Coulomb scattering explicitly included. We find 
that a finite steady-state terahertz spin polarization induced by the terahertz laser field, first pre- 
dicted by Cheng and Wu [Appl. Phys. Lett. 86, 032107 (2005)] in the absence of dissipation, exists 
even in the presence of all the scattering. We further discuss the effects of the terahertz laser fields 
on the spin relaxation and the steady-state spin polarization. It is found that the terahertz laser 
fields can strongly affect the spin relaxation via hot-electron effect and the terahertz-field-induced 
effective magnetic field in the presence of spin-orbit coupling. The two effects compete with each 
other, giving rise to non-monotonic dependence of the spin relaxation time as well as the amplitude 
of the steady state spin polarization on the terahertz field strength and frequency. The terahertz 
field dependences of these quantities are investigated for various impurity densities, lattice temper- 
atures, and strengths of the spin-orbit coupling. Finally, the importance of the electron-electron 
Coulomb scattering on spin kinetics is also addressed. 

PACS numbers: 72.25.Fe, 72.25.Rb, 71.70.Ej, 72.20.Ht 



I. INTRODUCTION 

Generating and manipulating spin coherence of elec- 
trons is one of the most important research focuses of 
semiconductor spintronics communityji'^'^ There have 
been many proposals to use electric field rather than 

magnetic field to generate and manipulate electron spin 
coherence4^i&i^>iaaiii2a3^a5i6^i8^ The mecha- 
nism of such proposals is that when the spin degree of 
freedom is coupled to the orbital degree of freedom via 
spatial varying (7-tensor (or magnetic field) or spin-orbit 
coupling (SOC) (such as the Rashba,^- the Dresselhaus^^ 
and the strain- induced^^ SOC), the electric field can act 
directly on spin through driving the orbital motion. Re- 
cently, Kato et al. achieved coherent spin rotation via 
gigahertz electric field applied along the growth direction 
of the g-tensor engineered GaAs/A^Gai-^^As parabolic 
quantum wellsi^ It has also been demonstrated exper- 
imentally that the SOC can enable electrical control of 
spin coherence without magnetic field i^'^'^^i^^ Rashba and 
Efros showed that even preferable to ac magnetic fields, 
ac electric fields can efficiently induce spin resonance in 
the presence of SOC in quantum wells, especially when 
the ac electric field is the in-plane one.^'^'*'^*' This effect 
is called electric dipole spin resonance which was later 
observed by Kato et al. in bulk GaAs," Meier et al. 
in GaAs/InGaAs quantum wells^ and Nowack et al. in 
GaAs quantum dots3i In these investigations, only weak 
electric fields are applied. 

Recently, Cheng and Wu showed theoretically that in 
InAs quantum wells, a strong in-plane terahertz (THz) 
electric field 1 kV/cm) can induce a large spin polar- 



ization (~ 10 %) oscillating at the same frequency of the 
THz driving field when dissipations are not consideredi^ 
This indicates that using strong THz electric field is a 
promising way to achieve high-frequency spin manipula- 
tion and spin generation in InAs-based nanostructures 
where the spin splitting is of the order of THzi^i^ 
They also showed that the strong THz field can greatly 
modify the density of states via the dynamical Franz- 
Keldysh effect,'^*' the sideband effect — and the ac Stark 
effect i^i^i^ii^ Later, Jiang et al. predicted similar ef- 
fects in singly charged InAs quantum dotsi^ As the den- 
sity of states of the electron spin system is greatly modi- 
fied by the intense THz fields, the dissipation effects may 
also be manipulated. However, up till now there is few 
study on the dissipative kinetics of strongly-driven elec- 
tron spin system, especially from a fully microscopic ap- 
proach. Previously, we have demonstrated that intense 
THz driving field in InAs quantum dots can elongate 
spin relaxation time (SRT) by more than one order of 
magnitudcj?^ The underlying physics is that the sideband 
effect strongly modulates the phonon-induced spin-flip 
transition rates. The effects of intense THz fields on spin 
relaxation in two dimensional electron system (2DES) 
are still unknown. The spin relaxation mechanism in 
2DES is quite different from that in quantum dots. In 
the driving-field-frce limit, it is widely accepted that spin 
relaxation in 2DES is dominated by D'yakonov-Perel' 
(DP) mechanismJ^ Previously, spin relaxation and spin 
dephasing have been closely studied using the kinetic 
spin Bloch equation approach developed by Wu et alM- 
in intrinsic, n-type and p-type semiconductors, in both 
Markovian and non-Markovian limits, even in systems far 
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away from the equihbrium (under strong static electric 
field or with high spin polarization) ^^liSiilii^ii^iMdSde 
The theory, which includes all relevant scattering (such 
as the electron-impurity, electron-phonon and electron- 
electron Coulomb scattering) explicitly, agrees very well 
with experiments i^^'^^ Many predictions from the theory 
have recently been confirmed experimentally.^*-^^ 

In this work, we first extend the theory to study the 
spin kinetics under intense THz laser fields in InAs quan- 
tum wells. The kinetic spin Bloch equations are de- 
rived in the spirit of the Floquet-Markov approacbSS via 
nonequilibrium Green function methodj^ The Floquet- 
Markov approach combines the Floquet theory, which 
solves the time-dependent Schrodinger equation of the 
strongly-driven system non-perturbatively^^^ with the 
Born-Markov approximation which is widely used in the 
derivation of the equation of motion for the reduced den- 
sity matrix of the concerned system.— The theory is 
frequently applied in the study of dissipative dynamics 
of strongly-driven systems i^i^i^^ With the extended ki- 
netic spin Bloch equations, we are able to investigate 
the effect of the strong THz fields on spin kinetics. We 
show that the steady-state THz spin polarization in- 
duced by the THz laser field, first predicted by Cheng 
and Wu in the dissipation-free case,— still exists in the 
presence of the full dissipation. Moreover, we investi- 
gate how this spin polarization as well as the spin relax- 
ation are manipulated by the external THz laser fields 
under various conditions. The predicted spin dynamics 
can be readily confirmed by Faraday and/or Kerr rota- 
tion measurement o^^'^^ under intense THz irradiation. 

The paper is organized as follows: In Sec. II, we set up 
the model and establish the kinetic spin Bloch equations. 
In Sec. Ill, we present our numerical results. We conclude 
and discuss in Sec. IV. 



II. MODEL AND FORMALISM 
A. Model 

We consider a 2DES confined in InAs quantum well. 
The confinement along the z direction is so strong (well 
width a = 5 nm) that only the lowest subband is con- 
sidered. The THz electric field ErHzit) = Esin(f^i) is 
applied in the quantum well plane. Here O = 2nv is 
the angular velocity, with v being the frequency of the 
THz field. In experiments, the THz field can be provided 
by the free electron laseri^ In the Coulomb gauge, the 
vector potential is A(t) = Ecos(r2i)/r2 and the scalar 
potential is = 0. In InAs quantum wells, the dominant 
SOC is the Rashba SOC.^o As the Rashba SOC is rota- 
tional invariant, we take Ethz along x axis. The total 
Hamiltonian is then 



H — He + Hei + He 



H. 



ep 



H 



ph • 



(1) 



Here Hph — X)a q '^^^Q'^l q^^.Q represents the phonon 
Hamiltonian (We take h = 1 throughout the pa- 



per). Fe. = EkQ,^j^Q^(«9^)e"*^'^^4^Ck-q<x, Hee 

5Ekk'Q,aa'^Ql^(*9^)l^4<TCk'<T'Ck'-q^'Ck+q<T and H^p 

EA,Qk^A,Q/(%)(aA,Q 



''A.-qi^kcr 



)cLck- 



denote the 



electron-impurity, electron-electron and electron-phonon 
interactions, respectively. Here, wa.q is the phonon 
frequency, q-a.q (cko-) is the phonon (electron) annihi- 
lation operator with A being the phonon branch in- 
dex (cr denoting the electron spin index), Rj stands 
for the position of jth impurity, Q = (q, 9z) is the 
three-dimensional momentum whereas q and k's are the 
two-dimensional ones along the well plane. I{iqz) = 

(^)2(g»9.a _ l)/{j^^Q[(^)2 „ ^2]} jg ^j^g foj.^ ^^^^^^ 

of the lowest subband. The matrix elements Vq, Uq 
for electron-electron and electron-impurity interactions 
as well as Ma,q for the electron-longitudinal-optical- 
phonon (electron-LO-phonon) interaction can be found 
in Rcf. 42 and the ones for the electron-acoustic-phonon 
interaction can be found in Ref.lH. We apply the random 
phase approximation in the screening of the Coulomb 
potential4i 

The electron Hamiltonian can be written as 



He = J2 H^'^' {k,t)cijy.,> 

k(T<T 



(2) 



where 



Ho(k,t) = iii±^^i 

' ' 2m* 



(7x 



- dyik^ + eA{t)) 

= Ek + TB^a;^^ C0s(17i) -I- i?em(l + C0s(2rit)) 1 

+ aii{(7xky — dykx) — aRO-yeE cos{flt)/Q,. (3) 



Here Ck = 7i3 



eE 

2m* ' IJ^ — m*Q,-^ ""^^^ ^em 

noted that the last term manifests that the THz electric 



and Eem = iliw- is 



field acts as a THz magnetic field along the y axis 
Sgff = 2aReEcos{m)/ {\g\^iB^) , 



(4) 



where g is the electron factor. We will show later that 
this THz-field-induced effective magnetic field has many 
important effects on spin kinetics. The term proportional 
to Eem is responsible for the dynamical Franz-Keldysh 
effect 1^1^ This term does not contain any dynamic vari- 
able of the electron system and thus has no effect on 
the kinetics of the electron system. Usually, the largest 
time-periodic term is the term ^Ekx^cos{Q,t), where the 
sideband effect mainly comes from. Under an intense 
THz field, this term can be comparable to or larger than 
ek- It should be noted that this anisotropic term breaks 
down kx —f —kx symmetry of the Hamiltonian. We will 
show later that this asymmetry leads to nonzero value 
of the average of kx over the electron system when the 
momentum scattering is included. 

The Schrodinger equation for electron with momentum 
k is 



zat^'k(t) = ^o(k,i)*k(t). 



(5) 
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According to the Floquet theory,^^ the sohition to the 
above equation reads 



*k„(i) 



0l(2:)^kr,(i) 



with 7y = ± denoting the spin branch and 0i(z) be- 
ing the wavefunction of the lowest subband. ^kjjl^) — 
g-iyk„t u,'^26'"^*Xcr where j/kr; and wJ^J are the eigen- 
value and eigen-vector of the equation 



+ a^iky + iak^)v^''^^. (7) 



This equation is equivalent to Eq. (2) in Ref. 19|. For each 
k, the spinors {|Ckr;(i))} at any time t form a complete- 
orthogonal basis of the spin spacci^i^ii^ The time evo- 
lution operator for state k can be written as 



X e 



'i[Eemt+Eem sin(20t)/(2Sl)] 



(8) 



B. Kinetic spin Bloch equations 



The kinetic spin Bloch equations offer a fully micro- 
scopic way to study spin dynamics in semiconductors, 
even in system with large static electric field where the 
hot-electron effect is important ..^ iM The electric field de- 
pendence of spin dephasing time in such system was stud- 
ied first by Weng et al^ for high temperature case and 
then by Zhou et al^ for low temperature case. In these 
works, the electric field appears only in the driving term. 
However, in the case of a strong time-periodic field, stud- 
ies have shown that including this field only in the driving 
term is insufficient <^ The correct way is to evaluate the 
collision integral with wavefunctions which are the so- 
lutions of the time-dependent Schrodinger equation, i.e., 
the Floquet wavefunctions, instead of the eigen wavefunc- 
tions in the field-free limit Moreover, the Markovian 
approximation should be made with respect to the spec- 
trum determined by the Floquet wavefunctionsi^ These 
improvements constitute the Floquet-Markov theory i^i^ 
Generally, this theory works well when the driven system 
is in dynamically stable regime and the system-reservoir 
coupling can be treated perturbatively. Besides giving 



good results, this approach has the advantage of be- 
ing easy to handle, compared with the rather compli- 
cated path-integral approach^ which makes it a use- 
ful approach in the study of spin kinetics under strong 
time-periodic fields. In this work, we incorporate the 
Floquet-Markov approach in setting up the kinetic spin 
Bloch equations. By making the Markov approximation 
with respect to the spectrum determined by the Floquet 
states, we first establish the kinetic equations for the sin- 
gle particle density operator. We then use the Floquet 
states as basis functions to expand the kinetic equations 
and obtain the kinetic spin Bloch equations in the pres- 
ence of the strong THz field. A similar approach has been 
applied to study the spin relaxation in singly charged 
quantum dots under intense THz driving fields in our 
recent work."^^ 

The kinetic spin Bloch equations for the single particle 
density operator can be written as^ 



(9) 



where dtp]si{t)\coh and dtp\i{t)\scat are the coher- 
ent and scattering terms respectively. Pk(i) = 

Sr,i,,2 Tr{cL,2Ck,n/5^(0}lm>(772| represent the 2 x 2 
single-particle density operators, with p'^{t) and {I??)} de- 
noting the density operator of the electron system and a 
complete-orthogonal basis in spin space separately. The 
explicit form of the equations without the intense driv- 
ing field can be found in the work of Cheng and Wu.— 
The coherent terms, which describe the coherent pre- 
cession determined by the electron Hamiltonian Hf, and 
the Hartree-Fock contribution of the electron-electron 
Coulomb interaction, can be written as 

dtMt)\coh = -z[i?o(k,t),/5k(t)] -z[S''^(k,<),pkW]. 

(10) 

Here I]«^(k,t) = - Ek',,. ^k-k',,J/fe)Ppk' W is 
the Coulomb Hartree-Fock self-energy. The scat- 
tering terms are composed of terms due to the 
electron-impurity (i9tPk|ei); electron-phonon (9tPk|ep), 
and electron-electron (9tpk|ee) scattering respectively. In 
the interaction picture, or the "Floquet picture" , 

P^{t) = K ^(k,t,0)pk(t)C/o'(k,i,0), (11) 

and under the generalized Kadanoff-Baym Ansatz;^ 
these scattering terms read 



dtpi^'''\t) 



E 



{pl'''''^\t)pr''^'^'^'\t) - p<r'^--\t)pl'^^^^\t)) } + {, ^ ,f " 



(12) 
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7r|MA,k-k',g, 



- A^A!k-k',.Pk'^^"""'(opr'''"'''(^)) } + - 



>F(r;ir;2) ^ ,n <F{Tq3v') 



(13) 



k',k", n,n' r;i...?y7 



(rpivVi)/^ n\rp(n'){V2n3)rp(n-n')(ri4rir,)rp{ri6V7) g\ 
I k,k' V ' / k',k k",k" — k~f-k' k" — k-|~k'.k" V ' / 

<F(ri3ri') 
k 

p<^(''^''^)(opr*'"'''Hi)Pk'^i'k+'kki)p^-^^''^''^'(i)) } + - 



£k"~k+k' 



.a ) (p^,^^"^"^^ (i)pr^''^'''^ (OPk'^ir+'k^ (OPk'^^"^"^^ (0 



(14) 



In these equations, A'^^k-k'.g, = ^A,k-k',g, + ^(1 ± 1) 
stands for the phonon number, is the impurity density, 
Pk = i - Pk, Pk = Pk, and Ekr, = £k + ^k,?- 

^&k^'(i.o) = (ek',,(t)ia.2W> 

X g«[(ek'-ek)t+7E sin(nt)(fe^~fe^)] 
^ ^ 5'^^)(')i'72)git(riO+ek,^^-£k„2)^ (15) 
n 

with 

k',k ~ / y a ^n+ma- \^^) 

ma 

Here F^''^ = X^m "^n+m aJmilskx) with J„i(a;) standing 
for the mth order Bessel fonction. 

= ^y^;'K')l')2)git("Sl+£k'r,i-^M2)^ (17) 

n 

with 

-"■k'jk ~ / , '^m a '-'n+m a- \-^°) 

ma 

{rj <-> 7]'} stands for the same terms as in the previous 
{} but with the interchange r] rj' . The term of the 
electron-electron scattering is quite different from those 
of the electron-impurity and electron-phonon scattering, 
as the momentum conservation eliminates the term of 
g27E sin(nt)fe^ ^ The coherent term in the Floquet picture 
reads 

9tP^ii)\cok = * E ^k-k',,J/(«g.)|' 

X [^k,k'(i,0)p^,(t)^k',k(i,0), P^(i)].(19) 

At zero THz field, the sideband summations are omitted 
and the above equations go back to those in Ref. 43. In 



Appendix A, we use the electron-impurity scattering as 
an example to show how to derive the scattering terms 
in the Floquet-Markov limit. 

The above equations clearly show the sideband ef- 
fects, i.e., nn in the (5-functions. The extra energy, 
nfl, is provided by the THz field during each scatter- 
ing process. This makes transitions from the low-energy 
states (small k) to high-energy ones (large k) become 
possible, even through the elastic electron-impurity scat- 
tering. These processes are the sideband-modulated 
scattering processes. For example, the weight of the 
n-th sideband-modulated electron-impurity scattering, 
I '5'k" k''^''^ ' I , is approximately 6,,^,^^ | J„ (7^ (/c^; - fc^ ) ) | when 
'-fEkx^cos{^t) in the Hamiltonian [Eq. ^] is the main 
source of the sideband effect. This term is important 
when n is around ±Nm, with representing the inte- 
ger part of jsikx — k'^). In fact, the sideband-modulated 
scattering makes the electron distribution in the three en- 
ergy ranges around e^, ek±N„i^ tend to be more uniform 
according to Eq. p^ . This, together with the other two 
scatterings, lead to the thermalization of the electron sys- 
tem, i.e., the hot-electron effect. Consequently, the elec- 
tron temperature Tg becomes larger than the lattice tem- 
perature T. Previously, it has been found that the hot- 
electron effect has large influence on spin dephasing and 
spin relaxation under high static electric field4^'^ In this 
paper, we will show similar effects in the case with THz 
driving field. Finally, it is noted that, as the sideband 
effect mainly comes from the term of "^e^x 

ricos(f2i), the 

electron-impurity and electron-phonon scattering plays 
the leading role in transferring energy from the THz elec- 
tric field to the electron system. 

A pronounced feature of the kinetic equations is that 
all the scattering terms arc directly time-dependent. In 
our previous study on spin dynamics in quantum dots 
with strong THz field)^ due to the fact that the spin- 
flip electron-phonon scattering rates are much smaller 
than the Zeeman splitting and the THz frequency, one 
can use the rotating-wave-approximation (RWA) treat- 
ment of the scattering terms and consequently only the 
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time-independent terms are kepti^ Here, as the scat- 
tering rate (especially that due to the electron-electron 
Coulomb scattering) is of the same order of the THz fre- 
quency, the RWA is no longer applicable. Thus the scat- 
tering terms become explicitly time-dependent. More- 
over, the scattering and coherent terms are time-periodic 
functions with period Tq = 27r/i7. Consequently the ki- 
netic spin Bloch equations are time-periodic differential 
equations, whose eigen-modes have the general form of 
~a ^ gi^£t Q^."e""* (a = 1, 2, 3, 4) according to Flo- 
quet theoremi^ Therefore, the solutions of the equations 
can be expressed as = J2a ^aP^ with Cq denoting the 
time-independent coefficients. 

The kinetic spin Bloch equations are solved numeri- 
cally with the numerical scheme laid out in Appendix B. 

After that, ^'^^ each k is obtained. From 

pf'^'^w = {Ci.,muo t(k,t,o)pk(t)?7„^(k,t,o)iek,'(o)) 

= (a,(i)l/5k(i)IWW), (20) 

by performing an unitary transformation, one comes to 
the single particle density matrix pk{t) in the coUinear 
basis {|(t)} which is composed by the eigen-states of a^. 
In this spin space, the spin polarization along any direc- 
tion can be obtained readily, e.g., Sz = J2k liPk ~ Pk)' 
= EkRelpJ,-^}, Sy - -EkMpJc-^}- From the tem- 
poral evolution of Sz, the SRT is extracted. 

Finally we briefly comment on the gauge invariance. 
Although the above formalism is derived in the Coulomb 
gauge, the obtained physical observables, e.g., Sz, is 
gauge invariant. This is because Sz = Tr(i(T2pk) = 
J2kaa' 5(kcr|o-2|kcr')(kcr'|/5k|kcr). Any gauge transforma- 
tion, A — !■ A -f Vx(i", and $ ^ $ — dtx{Y,t), gives 
|ka) e"*'^^^'"'*'' jkcr). Thus the results are gauge invari- 
ant. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

We numerically solve the kinetic spin Bloch equations, 
Eq. with all the scattering mechanisms explicitly 
included, to study the spin kinetics in InAs quantum 
wells under intense THz laser fields. The parameters 
used are listed in Table Ii^ The density of the 2DES is 
= 10^^ cm~^ and the quantum well width is a=5 nm 
throughout the paper. The Rashba parameter and the 
frequency of the THz field is taken to be afl=30 meV-nm 
and 1/ = 0.65 THz respectively, unless otherwise specified. 
The initial distribution of the electron system is chosen to 
be a thermalized distribution under the THz field, which 
is obtained by sufRcient long-time (typically ~10 ps) evo- 
lution from a spin-polarized Fermi distribution at the lat- 
tice temperature T: pl^ = 0, p^" = 1/ [e(^'=-^<')/'=«^ + l] 
(/icr denotes the chemical potential of electrons with spin 
tj) with the SOC being turned off.^^ 

The following of this section is divided into two parts. 
In the first part, we study the spin pumping due to the 



THz laser field. We first show that the THz field can 
pump spin polarization, first predicted by Cheng and 
Wu in the dissipation-free casc;^ even in the presence 
of full dissipation. We then investigate the amplitude of 
the steady-state spin polarization as function of the THz 
field strength and frequency for various impurity densi- 
ties, lattice temperatures and Rashba SOC parameters. 
In the second part, we investigate the spin dynamics with 
finite initial spin polarization. We first show the tempo- 
ral evolution of spin polarization for a typical case with 
different THz fields. We then study the dependence of 
the SRT on the strength and frequency of THz field under 
various conditions. 

TABLE I: Parameters used in the calculation 

D 5.9 X 10^ kg/m^ Vst 1.83 x 10^ m/s 

Vsi 4.28 X 10^ m/s ei4 0.35 x 10^ V/m 

S 5.8 eV UJLO 27.0 meV 

Ko 15.15 Koo 12.25 

g -14.7 m* 0.0239mo 



A. Spin pumping 

1. Temporal evolution of spin signals 

In Fig. [ll(a), we plot the spin polarization along the y 
axis, Sy, as a function of time when the initial spin polar- 
ization is zero for E = 0.5 (solid curve) and 1.0 (dotted 
curve) kV/cm. The spin polarizations along z- and x- 
axes are always zero. We also plot the THz-field-induced 
effective magnetic field, B^q [Eq. ([4])], as dashed curve in 
the figure. It is noted that the THz field pumps a large 
(several percent) spin polarization which oscillates at the 
same frequency with the THz field. This feature coin- 
cides with what predicted in the previous work where no 
dissipations are considered 3 Nevertheless, it is interest- 
ing to see that there is a delay of Sy with respect to the 
THz-field-induced effective magnetic field -B^g-, which is 
different from the dissipation-free case. The time depen- 
dence of 5*^ falls into the general form 

^.w = E^°"co«[""(*-*s)]- (21) 

n>0 

with and denoting the amplitude and the de- 
lay time respectively. The delay times are due to the 
retarded response of the spin polarization to the spin 
pumping caused by the THz field. This can be revealed 
by the following simplified analysis. Approximately, Sy 
satisfies the following equation, 

dtSy = - {Sy - Sy{t)) /t,, (22) 

where Sy(t) is the instantaneous equilibrium spin polar- 
ization induced by -Bgg(t) due to Pauli spin paramag- 
netism. The factor I/t^ represents the spin relaxation. 
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Under the initial condition Sy(0) 
the following solution 



0, the equation has 



Sy{t) 



(23) 



When the THz field is strong, the instantaneous equilib- 
rium spin polarization has the form of multi-frequency 
dependence: Sy{t) — J2n ^y^"''^'^ ^ demonstrated in 
Ref. As the effective magnetic field is in the form 
of cosine function, Sy{t) should be in the form Sy{t) = 
^„^Q 25'^ cos(nr2t), where Sy is real. The solution of 
Sy{t) at Ts is hence given by 



2S- 



^0 ^(nl7r,)2 + 1 



cos[nn{t ~ tn)], (24) 



with 



in 



arctan(?T,f2rs)/ (nfi). 



(25) 



Comparing the above equation with Eq. (HU, one ob- 
tains 5'°" = 25^/v^(nrJrl)2 + 1 and t"^ = tn- The delay 
time iJJ is indeed due to the retarded response of the spin 
polarization to the spin pumping. In the limit of Tg 0, 
one has Sy{t) = Sy{t), i.e., the spin polarization com- 
pletely follows the spin pumping due to the THz field. 
This is exactly the property of the results obtained in 
the previous dissipation-free studies j^ i^^'^*^ 

Typically n runs in the range of n = 1,2 in the pa- 
rameter regime of our investigation. For small THz field 
strength, E < 0.4 kV/cm, only the term with n — 1 con- 
tributes to Sy. For larger field strength, the term with 
n > 2 also contributes and the peak of Sy{t) is not sym- 
metric any more. Although the term with n > 2 may also 
contribute, the most important contribution still comes 
from n = 1 term. Consequently Sy signal still has good 
periodic behavior. 

It should be mentioned that under the RWA, the ki- 
netic spin Bloch equations in the Floquet picture is 
explicitly time-independent i ^6 , so Thus, the steady-state 
density operator in the Floquet picture becomes time- 
independent^*^ and the time-dependence of the spin po- 
larization Sy in the steady state becomes totally deter- 
mined by the time-evolution of the Floquet states, which 
completely follow the spin pumping due to the THz field. 
As a result, the RWA loses the important information of 
the retardation of the spin polarization to the THz field. 
In our study, we go beyond the RWA. 

As pointed out after Eq. ^ that the anisotropic term 
7£;fca;f2 cos(f2i) in the Hamiltonian breaks down — > 
—kx symmetry. However, without scattering, the den- 
sity matrix in the Floquet picture does not change with 
time and is determined solely by its initial value. For 
the choice of isotropic initial distribution, the k^ — > —k^ 
asymmetry of the density matrix never shows up. In the 
presence of scattering, the density matrix should show 
the asymmetry of the Hamiltonian. Consequently, the 
average of k^ of the electron system, (k^), is nonzero. 
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FIG. 1: (Color online) (a) Spin polarization along y axis, Sy, 
as function of time with zero initial spin polarization for E = 
0.5 kV/cm (solid curve) and 1.0 kV/cm (dotted curve). T = 
50 K and Ni = 0.05Ne. The dashed curve is the THz- field- 
induced effective magnetic field -Bgff with i5 = 1.0 kV/cm. 
Note that the scale of the dashed curve is on the right hand 
side of the frame, (b) Sy vs. time for _E = 1.0 kV/cm with 
Bgjj included (solid curve) and excluded (dotted curve). 



Below we will show that, quite remarkably, the scattering 
terms within the RWA keeps the symmetry of k^ —k^ 
and the scattering terms which do not keep the symme- 
try only appear in the time-dependent (beyond RWA) 
scattering terms. 



Look at, e.g., the electron-impurity scattering 
[Eq. p^ ]. the weight of the n-th sideband-modulated 

scattering is P^""!"^"^) = 5^;^P(t, 0)4?^^"^"^'. As 
the main source of the sideband effect is the term 
7£;fca;il cos(f2i) in the Hamiltonian [Eq. the weight 
is approximately 



pinmmm) 



S-QTIi '^772 1)3 



'[(Ek— ek')*+7-B sin(nt)(fcx— fe^)] 



(26) 



Pn can be further decomposed into the time-independent 
P™ and the time-dependent = P^^ + P^"^ parts (omit- 
ting the superscripts of 77): 



pi 



Jn {1e{K " ^x)) Jn ilEiK ~ ^x)) 



(27) 
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m-.even m^o 

X J„[7£(fc;-A:,)]e-""^*, (28) 

Pf = E Jn+mbEiK-k^)] 
m:Odd 

X J„[7ij(fc;-A:.)]e-"""*. (29) 

It is seen that under the transformation: k'^ —> —k'^ and 
kx —kx, and are invariant but P^^ is changed. 
This indicates that a portion of the time-dependent (be- 
yond RWA) scattering terms break down the kx —kx 
symmetry. With these scattering terms, the density ma- 
trix should evolve to be asymmetric in kx direction. This 
leads to (kx) ^ 0. (k^) should also oscillate with time as 
P^^ does. 

In the presence of the SOC, (kx) leads to a second 
effective magnetic field: 

Bay{t)^2aB.{kx)/{\g\fiB)- (30) 

Indeed, we find that the spin polarization Sy is still 
nonzero when B^g is turned off by omitting the corre- 
sponding term in the Hamiltonian. In Fig.[T]^b), we plot 
Sy{t) for both cases with and without B^q. It is seen 
that Sy is nonzero when Pgff is excluded, although the 
amplitude is reduced. This spin polarization is induced 
by Bav via Pauli paramagnetism. The results indicate 
that Bav oscillates with time and is smaller than Pgff. 
Moreover, there is a change in the delay of the oscilla- 
tion due to different time-dependence of Bav compared 
to Pgff. This difference also contributes to the delay of 
the spin polarization which is induced by the total effec- 
tive magnetic field B {B = Pgg: -I- Pai>)- 

Finally, it is found that a small initial spin polarization 
(~ 4%) along the z axis makes marginal effect on the time 
dependence of Sy. 

2. Steady-state spin polarization 

In this subsection, we discuss the dependence of the 
amplitude of the steady-state spin polarization Sy (the 
peak value of Sy) on the THz field. 

In Fig. [5] we plot the amplitude of the steady-state 
spin polarization (ASSSP) as a function of THz field 
strength for the cases with and without the THz-field- 
induced effective magnetic field Pgff. Two typical lat- 
tice temperatures T = 50 K and 100 K are investigated 
with the impurity density Ni = Q. It is seen that the 
ASSSP first increases then decreases with the strength 
of the THz field. The effective magnetic field increases 
with the THz field strength. According to Pauli para- 
magnetism, however, the spin polarization should always 
increase with the magnetic field. Here the decrease of the 
ASSSP mainly originates from the hot-electron effect. To 
elucidate this point, we plot the hot-electron temperature 
Tf. in Fig. [3] (the method used to obtain is given in Ap- 
pendix C). It is seen that the hot-electron temperature 
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FIG. 2: (Color online) Dependence of the ASSSP on THz 
field strength for T = 50 K (solid curve with •) and 100 K 
(solid curve with □) without impurities {Ni — 0). The dotted 
curves are the same as the solid ones but without the THz- 
field-induced effective magnetic field Sgg-. 



increases with the THz field strength. The increase of the 
hot-electron temperature decreases the induced spin po- 
larization according to Pauli paramagnetism. It is noted 
from the figure that the largest ASSSP can be ~ 3.5x10^ 
cm~^ which corresponds to a large spin polarization of 
7 %. This indicates that the intense THz field is a very 
efficient tool in generating spin polarization. It can be 
noticed in Fig. [2] that the ASSSP is smaller at higher 
temperature. The decrease of the ASSSP is due to the 
increase of the hot-electron temperature with the lattice 
temperature, as indicated in Fig. [31 
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FIG. 3: (Color online) Hot-electron temperature Te as func- 
tion of THz field strength for T = 50 K (solid curve with •) 
and 100 K (dashed curve with □) without impurities {Ni — 0). 
The dotted curve with • is the same as the solid one, but with 
only n = 0, ±1 allowed in the sideband-modulated scattering. 
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It should be pointed out that, differing from our pre- 
vious study on spin relaxation in quantum dots,"^^ here 
the n-th sideband-modulated scattering rate differs little 
from each other. The energy conservation (the S func- 
tions in the scattering terms) gives different final state 
for different n with given initial state, thus the momen- 
tum transferred into the system can change effectively 
with n. However, the matrix elements of all the scat- 
tering mechanisms vary slowly with the momentum due 
to the screening and the quantum confinement along the 
growth direction. Consequently, the sideband-modulated 
scattering rate varies slowly with n and the manipulation 
of the spin relaxation via sideband modulation of the 
spin-flip scattering does not apply in 2DES. In 2DES, 
the main effect of the sideband-modulated scattering is 
the hot-electron effect. As we have pointed out, the n-th 
sideband-modulated scattering tends to make the distri- 
bution be flatter in the energy range of nQ, which thus 
leads to the hot-electron effect. In Fig. [31 we also plot 
the hot-electron temperature when the summations of 
n in the scattering terms [Eqs. ([T^ . (fT3)) and IH)) ] are 
restricted to n = 0,±1. In the previous studies on the 
effect of THz field on spin dynamics, only these processes 
are considered, where the THz field is weak j^^i^^'^" It is 
seen that the hot-electron temperature is largely reduced 
by the number of sideband involved in the scattering, 
especially when the THz field is strong and hence the 
sideband- modulated scattering with \n\ > 1 is important. 
This confirms the important role of sideband-modulated 
scattering to the hot-electron effect. 
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FIG. 4: (Color online) Dependence of the ASSSP on THz 
field frequency for E — 1.5 kV/cm (•) and 2.5 kV/cm (□). 
T = 100 K and A^^ = 0. 

The dotted curves in Fig. [21 is the ASSSP calculated 
without the THz-field-induced effective magnetic field 
Beff. As analyzed before, here the ASSSP is induced by 
Bav The contribution of Bav becomes more important 
when THz field strength increases. This is because that 
the terms P,f^ increase with the sideband effect which 



increases with THz field strength. Moreover, it is seen 
that the ASSSP due to Bav is larger at higher temper- 
ature (100 K in the figure) when the THz field strength 
is small. This is because the scattering terms leading 
to the breakdown of — > —kx symmetry increase as 
the electron-LO-phonon scattering is more efficient at 
T — 100 K. However, when the THz field strength is 
larger, the hot-electron effect becomes more important. 
(We find that the hot-electron temperature changes lit- 
tle when BgQ- is removed.) The hot-electron effect also 
reduces the ASSSP induced by Bav At large THz field 
strength this effect becomes more important and the dif- 
ference of ASSSPs at 50 and 100 K becomes marginal. 
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FIG. 5: (Color online) Hot-electron temperature Te as func- 
tion of the THz frequency for E = 1.5 kV/cm (•) and 2.5 
kV/cm (□). r = 100 K and Ni = 0. 

We now turn to investigate the dependence of the 
ASSSP on THz frequency. It is noted that the am- 
plitude of the THz-field-induced effective magnetic field 
B^Q [Eq. ^] decreases with THz frequency. However, 
the hot-electron effect induced by the sideband effect 
which increases with je^ = eE / {m*^l), also decreases 
with THz frequency. These two effects again compete 
with each other. In Fig. [4] we plot the ASSSP as func- 
tion of the THz frequency for two cases: E — 1.5 and 
2.5 kV/cm with T = 100 K. It is seen that for the case 
with E = 2.5 kV/cm, the ASSSP first increases then de- 
creases with the THz frequency due to the competition 
of the two effects. To examine the hot-electron effect, 
we also plot the hot-electron temperature in Fig. [51 It 
is seen that the hot-electron temperature decreases with 
the THz frequency. For large THz frequency, the hot- 
electron effect is marginal and thus the ASSSP decreases 
with the THz frequency as the THz field induced effec- 
tive magnetic field does. For smaller THz frequency the 
hot-electron effect becomes dominant and the ASSSP in- 
creases with the THz frequency. Consequently there is a 
peak frequency where the ASSSP reaches the maximum. 
For the case with E = 1.5 kV/cm, there should be a peak 
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with the peak frequency being much smaller than the fre- 
quency we calculated. This is because the hot-electron 
effect is much weaker than the case with E = 2.5 kV/cm 
(see Fig. [S]). 



B. Spin dynamics with finite initial spin 
polarization 

1. Temporal evolution of the spin signals 

In Fig. [SI we plot the temporal evolutions of the spin 
signals along z, x and y axis at different THz field 
strengths. The initial spin polarization is taken to be 
4 % along the z axis. In Fig. EJa), one finds that Sz 
exhibits oscillatory decay. This resembles the low tem- 
perature spin decay observed in Refs. and [6ll, which 
is due to the large spin-orbit effective magnetic field and 
weak scattering)^ i.e., the system is in the weak scatter- 
ing limit. It is noted that Sz decays faster when the THz 
field strength increases. Moreover, the spin oscillation 
frequency also increases as indicated by the left-shift of 
the peak around 0.7 ps, which is due to the total effec- 
tive magnetic field B induced by the THz field. From Fig. 
[6jb) , one notices that a small value of Sx is excited but 
eventually decays to zero. This is again due to the effec- 
tive magnetic field B, which rotates Sz to Sx- The first 
peak value of Sx increases with the THz field strength. 
Without THz field, Sx = 0. In Fig. ^c), it is seen that 
Sy is also induced and reaches a non- vanishing oscillatory 
value after ~ 3 ps evolution, similar to what observed in 
Fig. [T] where the initial spin polarization is zero. 
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FIG. 6: (Color online) Temporal evolution of the spin signals 
(a): Sz; (b): Sx', and (c): Sy for THz field strengthes: 
kV/cm (solid curves), 1 kV/cm (dashed curves), and 2 kV/cm 
(dotted curves). Ni = 0.05Ne and T = 50 K. 



2. SRT 

In Fig. [71 the SRT, which is extracted via fitting the 
exponential decay of the envelop of Sz , is plotted as func- 
tion of THz field strength for different impurity densities 
Ni — 0, 0.02iVe, and 0.05Ne, with lattice temperature 
T = 50 K. 

We first discuss the case with Ni=0 (solid curve with 
• ). It is noted that the SRT first increases then decreases 
with the THz field strength. The underlying physics is 
that there are two consequences of the THz field: (i) the 
total THz field-induced effective magnetic field B; (ii) 
the hot-electron effect. Effect (i) can give a magnetic 
field as large as several Tesla (2.6 T per I kV/cm THz 
field with v — 0.65 THz, of which the corresponding Zee- 
man splitting is as large as 2.2 meV). This effective mag- 
netic field blocks the inhomogeneous broadening from the 
Rashba SOC. It thus elongates the SRT.^>29 rj.^^ ^^^^ 
consequences of effect (ii) are the enhancement of mo- 
mentum scattering as well as the inhomogeneous broad- 
ening as the electrons distribute on larger k states where 
the SOC is larger. Enhancement of the inhomogeneous 
broadening shortens the SRT according to our previous 
studieSfS Et ^ i ^ It is found that in the strong scattering 



limit, SRT increases with the momentum scattering.— i2i 
However, it is demonstrated in Rcf. 40 that in the weak 
or intermediate scattering limit, the SRT decreases with 
increasing the momentum scattering. In our case, due to 
the large Rashba SOC parameter, the system is in the 
weak/intermediate scattering regime. This can be fur- 
ther checked by the fact that the SRT decreases when the 
electron-impurity scattering is strengthened by increas- 
ing the impurity density, as shown in Figs.[7|and[8l Thus, 
effect (ii) shortens the SRT. The two effects compete with 
each other and hence the SRT varies non-monotonically 
with the THz field strength: For small THz field strength 
the increase of the total THz field-induced effective mag- 
netic field B is dominant. As a result, the SRT increases; 
For large field strength, the hot-electron effect becomes 
more important. Consequently the SRT decreases. 

To further elucidate the influence of effect (i), we re- 
move part of the THz field-induced effective magnetic 
field, B^Q, by excluding the term aj^aycE cos{rit)/fl 
from the Hamiltonian, and then calculate the SRT. We 
plot the obtained SRT as dotted curve in Fig. [71 It is 
seen that the SRT is reduced, especially at large THz 
field strength. It is checked that the hot-electron effect 
changes little when B^q is removed as it is not the main 
source of the hot-electron effect. The results confirm that 
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FIG. 7: (Color online) Dependence of the SRT r on THz field 
strength for impurity densities: iVi = (•); iV; = 0.02iVe (□); 
Ni = 0.05Ne (A). Solid curves: from full calculation; Dotted 
curve: from the calculation without the THz-field-induced 
effective magnetic field -Bgfj. 



0.7 




0.2 ' ■ 

12 3 



E (kV/cm) 

FIG. 8: (Color online) Dependence of the SRT r on THz field 
strength for impurity densities: Ni = (•); Ni = 0.02Ne (□); 
and Ni = 0.05iVe (A). T = 100 K (solid curves) and 50 K 
(dashed curves). 
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FIG. 9: (Color online) Dependence of the hot-electron tem- 
perature Te on THz field strength for different impurity den- 
sities: Ni = {• ); Ni = 0.02N^ (□); N, = 0.05iVe (A). 
T = 100 K (sohd curves) and T = 50 K (dashed curves). 
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FIG. 10: (Color online) Dependence of the SRT r on THz 
frequency for impurity densities: Ni = (•); Ni = 0.02A'^e 
(□); and Ni = 0.05A^e (A). Solid curves: from full calcula- 
tion; Dotted curve: from the calculation without the THz- 
field-induced effective magnetic field Bgff. E = 1 kV/cm and 
T = 50 K. 



the THz field-induced effective magnetic field indeed in- 
creases the SRT. 

For the cases with N, = 0.02Ne and 0.05iVe, the SRT 
decreases with the THz field strength monotonically. It 
is noted in Fig. 9 (solid curves) that the hot-electron 
temperature is larger at higher impurity density. The 
enhancement of the hot-electron effect overcomes the in- 
crease of the effect of the THz field- induced effective mag- 
netic field in these two cases, which leads to the mono- 
tonic decrease of the SRT. 

It should be mentioned that in the case of static electric 
field, the hot-electron effect is more important at smaller 



impurity density under a given electric fieldi i '^^ How- 
ever, under intense THz field the hot-electron effect is 
more pronounced at larger impurity density where the 
sideband-modulated scattering, which can transfer the 
THz photon energy into electron system, is stronger. 
Similar effects have been reported by Lei in the study 
of charge transport under intense THz fieldj^ 

We further discuss the temperature dependence of the 
SRT. In Fig. i] we plot the SRT as function of THz field 
strength at T = 100 K and 50 K. One can see that the 
SRT at high temperature (T = 100 K) is much smaller 
than that at low temperature {T = 50 K). This is because 
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FIG. 11: (Color online) Dependence of the hot-electron tem- 
perature Te on THz frequency for impurity densities: Ni — 
(•); Ni = 0.02iVe (□); and Ni = 0.05N^ (A). E = 1 kV/cm 
and T = 50 K. 



FIG. 12: (Color online) Dependence of the SRT r on THz field 
strength for different SOC parameters: an — 1 meV-nm (•); 
Or = 10 meV-nm (□); and ur = 30 meV-nm (A). T = 50 K 
and Ni = 0.05A^e. Note that the scale for the curves with □ 
and A are on the right hand side of the frame. 



the electron-LO-phonon scattering at high temperature 
is much more efficient than that at low temperature. The 
increase of scattering thus enhances the hot-electron ef- 
fect (see Fig. 9) and the enhancement of the hot-electron 
effect reduces the SRT. It is also noted that for all three 
impurity densities at 100 K, the SRTs decrease with the 
THz field strength monotonically. This indicates that 
the hot-electron effect of the THz field is dominant as 
the momentum scattering is strong. 

We now turn to the THz frequency dependence of the 
SRT. As has been demonstrated before, both the THz 
field-induced effective magnetic field and the hot-electron 
effect decrease with the increase of the THz frequency 
and they compete with each other on spin relaxation. 
In Fig. [TOl we plot the SRT as function of the THz fre- 
quency for three different impurity densities: A'^; = 0, 
0.02iVe, and 0.05iVe. The lattice temperature is T = 50 
K. It is noted that for the impurity-free case, the SRT 
first increases then decreases with the THz frequency, 
which is similar to the dependence of the SRT on the 
THz field strength. This is again due to the competi- 
tion between the hot-electron effect and the THz field- 
induced effective magnetic field. Similarly, for the cases 
with Ni = 0.027Ve and 0.05iVe, the SRT increases mono- 
tonically with the THz frequency. When B^q is removed, 
the SRT for impurity-free case is reduced and the peak 
frequency where the SRT gets maximum becomes larger. 
This indicates the weakening of the THz field-induced ef- 
fective magnetic field since the hot-clcctron effect changes 
little. 

Finally, we discuss the dependence of the SRT on the 
Rashba SOC parameter. In Fig. [T^ we plot the SRT as 
function of THz field strength at different Rashba pa- 
rameters, an = 1, 10, and 30 meV-nm. The lattice tem- 
perature is taken to be T = 50 K and the impurity den- 
sity is Ni = Q.05Ne. It is seen that for small Rashba 



SOC coefficient an = 1 meV-nm, the SRT first increases 
then decreases with the THz field strength. However, for 
large SOC, the SRT decreases monotonically with the 
THz field strength. As has been revealed previously that 
in the presence of the impurity density Ni = 0.05A''e, the 
hot-electron effect dominates the SRT. The hot-electron 
effect leads to the increase of both the scattering and the 
inhomogeneous broadening. For the case with — 1 
mcV-nm, which is in the strong scattering regime, in- 
crease of scattering leads to longer SRT; whereas the 
increase of inhomogeneous broadening leads to shorter 
SRT. Therefore, the two effects compete with each other: 
the SRT first increases due to the enhancement of scat- 
tering and then decreases due to the increase of inho- 
mogeneous broadening. This behavior is similar to the 
case under a strong static electric field4^ For the cases 
with larger SOC, which is in the intermediate scattering 
regime, both effects decrease the SRT. 



C. Effect of electron-electron Coulomb scattering 

Previously it has been shown that the electron-electron 
Coulomb scattering plays an important role in spin re- 
laxation due to the DP mechanism ^^iiSii^ii^iiiiHies p^j.- 

thermore, for systems under strong electric field, the 
electron-electron Coulomb scattering is crucial for the 
electron system to establish its quasi-equilibrium state.— 
Here we demonstrate that it also has nontrivial effects on 
the spin dynamics under intense THz field. In Fig. [121 
we plot the temporal evolution of the spin signals, \Sz\ 
and Sy, calculated with the electron-electron scattering 
included (solid curve) and excluded (dashed curve) un- 
der same initial distributions and conditions. It is seen 



12 















fYl 






w n 
» # ■ 
1 1 


(a) 











01234567 



t (ps) 



0.02 
0.01 

i 

3 
o 

>. -0.01 
1/1 

-0.02 











/ \ y\\ 


''M if V 




r \ 
1 


i'll J/A J^.l ,f\l 'Til .TM ! 

'/', f'. '/ ' ' / ' ', ' H ' 


r/ V'] 


V'l 




■ 


;l , 1/ 1 ' ' ' '1 1 1 ' 'I 

1 'I M i ' '\ ' '\ 1 i\ M M ' 
i 1 V'l ri VI V'l vl 'VI 


(b) 











5 10 15 20 



t (ps) 

FIG. 13: (Color online) Temporal evolution of the spin signals 
calculated with the electron-electron scattering is included 
(solid curve) and excluded (dashed curve), (a): \Sz\; (b): Sy. 
E = 1 kV/cm, T = 50 K and N^ = 0. 

from Fig. [131 a) that in the case with the electron-electron 
Coulomb scattering, \Sz\ exhibits good exponential de- 
cay, superimposed by the THz oscillations. Otherwise, 
the decay is non-exponential and the decay rate becomes 
much slower, which indicates that the spin relaxation is 
markedly reduced. It is further seen from Fig. fTSF b) that 
with electron-electron Coulomb scattering, Sy reaches 
the steady state much faster with a larger peak value. 
All these demonstrate the importance of the electron- 
electron Coulomb scattering to the spin dynamics. 



IV. CONCLUSION AND DISCUSSION 
A. Conclusion 

In conclusion, we have developed the kinetic spin Bloch 
equations for 2DES with Rashba SOC under intense 
THz laser fields, with all the relevant scattering mecha- 
nisms such as the electron-impurity, electron-phonon and 
electron-electron Coulomb scattering explicitly included. 
The formalism is very general and can be applied to study 
spin kinetics in many-body electron or hole system under 
strong time-periodic driving fields with arbitrary SOC. 
Moreover, our formalism goes beyond the RWA treat- 
ment of the scattering. By solving the kinetic spin Bloch 
equations numerically, we investigate the effect of the in- 
tense THz fields on the spin kinetics. We focus on the 
THz field-induced steady-state spin polarization and the 
effect of the THz field on spin relaxation. 

We first study the temporal evolution of the spin po- 



larization under intense THz field at zero initial spin po- 
larization. We find that the THz field can pump a fi- 
nite steady-state THz spin polarization in the presence of 
all relevant scattering. The spin polarization is induced 
by the THz field-induced effective magnetic field in the 
presence of SOC. The maximum spin polarization in the 
steady state can be as large as 7 %, which shows that 
the intense THz field is a very efficient tool in generating 
spin polarization. 

As our approach goes beyond the RWA treatment of 
the scattering, we find some interesting features which 
are absent in the RWA treatment. The first feature is 
that there is always a retardation of the spin polarization 
in response to the THz field-induced effective magnetic 
field. Another feature is that, as the Hamiltonian breaks 
the kx — *■ —kx symmetry via the term "fEkx^cos(flt), 
the average of kx over the electron system (kx) becomes 
nonzero and oscillates with time. In the presence of SOC, 
{kx) leads to another effective magnetic field which also 
induces spin polarization. We find that, remarkably, un- 
der the RWA, the kx —kx symmetry is still kept. 

We further study the dependence of the amplitude 
of the steady-state spin polarization on the THz field 
for different lattice temperatures, impurity densities and 
Rashba SOC parameters. It is found that the main con- 
sequences of the THz field are: (i) the hot-electron effect 
due to sideband-modulated scattering and (ii) the THz 
field-induced effective magnetic field due to the SOC. 
Both effects increase with the THz field strength but de- 
crease with the THz frequency. The amplitude of the 
steady-state spin polarization increases with effect (ii), 
but decreases with effect (i) according to the Pauli para- 
magnetism. At small THz field strength (and/or large 
THz frequency) the hot-electron effect is weak and ef- 
fect (ii) dominates. The amplitude of the steady-state 
spin polarization thus increases (decreases) with the field 
strength (frequency). At large THz field strength (low 
THz frequency), the hot-electron effect becomes more im- 
portant than effect (ii) and the amplitude of the steady- 
state spin polarization decreases (increases) with the THz 
field strength (frequency). 

We also find that the THz field can strongly change 
the SRT due to the two effects addressed above. Specif- 
ically, the hot-electron effect shortens the SRT via en- 
hancement of momentum scattering and inhomogeneous 
broadening for the system in weak/intermediate scatter- 
ing limit. Meanwhile, effect (ii) increases the SRT due 
to the blocking of the inhomogeneous broadening. For 
small impurity densities at low temperature, when the 
THz field strength is small (and/or the THz frequency 
is large), the hot-electron effect is weak, and effect (ii) 
becomes dominant. The SRT thus increases (decreases) 
with the THz field strength (frequency). At large THz 
field strength (small THz frequency) the hot-electron ef- 
fect is more important than effect (ii). The SRT thus 
decreases (increases) with the THz field strength (fre- 
quency). However, for large impurity densities or high 
temperatures, the enhancement of the hot-electron effect 
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overcomes the increase of the effect (ii). Consequently 
the SRT decreases (increases) with the THz field strength 
(frequency). We also discuss the SOC dependence of the 
SRT at large impurity densities where the hot-electron 
effect dominates. For small SOC, which is in the strong 
scattering regime, the SRT first increases with the THz 
field strength due to enhancement of momentum scat- 
tering, then decreases with it due to enhancement of 
the inhomogeneous broadening. For large SOC, which is 
in the weak/intermediate scattering regime, increase of 
scattering also reduces the SRT. Consequently, the SRT 
decreases monotonically with the THz field strength. 

B. Discussion 
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Finally we compare our study with the electric dipole 
spin resonance (EDSR) in the literature. To simplify 
the discussion, we introduce a simple spin Hamiltonian 
which characterize the spin dynamics of our Hamiltonian 
[Eq. ©] and the EDSR: 

Hsptn{t) = i(Ao + Ai)(Tj, + iA2(Tj^ -l-7Cos(m)(Tj,. (31) 

Here Ai and A2 characterize the k-dependent effec- 
tive magnetic fields due to the SOC. Aq represents the 
external static magnetic field used in the EDSR set- 
^p^i^^jA (For our case: Aq = 0, Ai = 2aRky, 
A2 = -2aRkx, 7 = -aneE/n.) In EDSR, Aq is usu- 
ally much larger than Ai, A2 and 7, and Ao — ^l^^^ 
To the lowest order approximation, the spin dynam- 
ics is governed by i?"pj„ = ^Aoax + j cos{D,t)(7y. In 
the RWA the solutions of the Schodinger equation are 
given by = e^'^*^{e-'^*/^x+ ± ie'^'/^x-) with 
axX± = ix±- With initial condition ■0 = X+: one ob- 
tains {Sx) = |cos(7i), which is the well-known Rabi os- 
cillation. The k-dependent Ai and A2 effective magnetic 
fields lead to the damping of the Rabi oscillation due to 
the DP mechanism in the presence of scatteringjii For 
the case of strong driving field (7 > i^), the solutions of 
the Schodinger equation are the Floquet wave functions 
^'kjj given in Eq. ([6]) . Now the spin dynamics is given by 

X g«*fek,-yk„'+("'-")f^]^ (32) 

with i = x,y,z. From the above equation, it is seen 
that, unlike the weak driving-field case where only a sin- 
gle Rabi frequency is observable, here the spin signal Si(t) 
oscillates at many frequencies ykn — ykri' + {m — n)fl (with 
771 — n = 0, ±1, ±2, • • • ). Moreover, in our case, Ai, A2, 7 
and f2 are on the same order of magnitude while Aq = 0. 
Thus the spin precession frequency varies largely with 
k. This large inhomogeneous broadening of spin preces- 
sion frequency smears out the driving field-induced Rabi 
oscillation of the spin polarization signals. 



APPENDIX A: DERIVATION OF 
ELECTRON-IMPURITY SCATTERING TERM 
IN FLOQUET-MARKOV LIMIT 

Here we give terms due to electron-impurity scattering 
as an example. Terms due to other scattering can be 
obtained similarly. From nonequilibrium Green function 
theory^ the electron-impurity scattering term can be 
written as 

dtPi^L = { A^^{><) + A^{<>)} + {...}', (AI) 
where 

X / dTU^ik',t,T)p>iT)p<iT)U§ik,T,t). iA2) 

A^{<>) can be obtained by interchanging > and <. It is 
better to work in the interaction picture, or the "Floquet 
picture" : 

P^W = t(k,<,0)p^(t)J7o^(k,i,0). (A3) 
After this transformation, the term becomes 

A^{><) = u;;\k,t,o)At\><mk,t,o) 

= E "'^k-lCgJ^l^?.)!' / dTSuM'it,0) 

xp>^{r)S^,^u{T,0)p<^{r), (A4) 

where S'k,k'(i,0) = Uq \k,t,0)U^{k' ,t,0). According 
to Floquet-Markov theory, the Markov approximation 
should be made with respect to the spectrum determined 
by the Floquet wavefunctions, i.e., p^ (r) « p^ {t). 
Thus, the scattering term becomes 

A^i><) = E"*'^k-k',,J^fe)l'^k,k'(i,0) 

k',g. 

xp^,^(i) /* dr5k',k(r,0)p<^(0. (A5) 
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The next step toward the expUcit form of the scattering 
term is based on the analysis of the elements of S. Ex- 
panding the kinetic equations in the basis of {|^kr;(0))}, 
the elements of S are given by 



X e 



»[(Sk' ^'^k)t+7B sin(nt)(fc^-fc^)] 



.(A6) 



The scattering term can then be explicitly laid out by 
expanding all the operators in the basis of {|^k))(0))}: 

k',(;j,n,77ir;2>)3 



X S(nfl + £k'r,2 ^ ^k 



(A7) 



According to Eq. (Al), one can readily arrive at the 
full expression of the electron-impurity scattering term, 
which is exactly Eq. p2)) . 



APPENDIX B: NUMERICAL SCHEME 

Our numerical scheme is based on the scheme laid out 
in detail in Ref. [i^ . where the nonlinear kinetic spin 
Bloch equations are solved self-consistently with high 
accuracy^i^i^i^ The scheme is based on a discretiza- 
tion of the two dimensional momentum space with Nx M 
control regions where the k-grid points are chosen to be 
k;^m = ^/2m*Ei (cos 9 rriT sin 9 rn)- In principle, the co- 
herent terms are easily solved. However, the scattering 
terms are difficult to solve as the (5-functions are hard 
to be integrated numerically. To facilitate the evalua- 
tion of the ^-functions in the scattering terms, we set 
El — {I + 1/2) AE and lulq — tllo^E where I and 
TiLO are integer numbers and AE is the energy span in 
each control region4^ To apply this scheme to the ki- 
netic spin Bloch equations with THz field, we also set 
ri = TT-THzAi? with tithz being integer number (typically 
1 ~ 3 in our calculation). However, the 5- functions 
are still difficult to be evaluated as j/k?; and AE are 
not commensurable. We therefore use the approxima- 
tion j/kr) ~ -^k??^-^' ^it'^ -^kr; b^iug the integer part of 
y\^ri/AE. This approximation affects the spin kinetics 
marginally as |yk»)— -^k^Ai?| is usually much smaller than 
ksT and/or the chemical potential. Moreover, as the 
driving field is very strong, the spectrum of the Floquet 
states is mainly determined by the sideband effect and 
2/kr; — N^^AE only plays a quite marginal role. Further- 
more, one can approach the exact results by increasing 
TT-THz- In our computation, we make sure that for n^Hz we 
choose, the relative error is less than 5 %. To make the 
treatment consistent, we also approximate S^/-^^\t, 0) « 



J2n ^k' 



{n){mn2)jtlnn+e^,-e^ + {Ny,^^-Ny^^)AE] 



Or more 



concisely, 5k',k(t,0) « E„ ^L" ke'*^"^^+"-'""' 

7^k^k(^,o), where 7^^7)^''^''^) = si", 



(m){r)iri2) 



k',k 



with m sat- 



isfying mriri 



Nl 
k'771 



y _ 

kri2 



Correspondingly, 



rk',k(i,0) « ^^wH^e^tCnAB+e^-ek) = W^,^^(t,Q), 

with W^T^J"^"^^ = yM(rM2)_ j^ggp ^j^g coherent pre- 
cession due to 2/k?) — N^T^^E by adding it into the coher- 
ent term. After these approximations, the coherent and 
scattering terms of the kinetic spin Bloch equations read 

<9*P^WL„;,-4 ^ ^k-k',gJ/fe)PWk.k'(i,0)p^,(i) 

X Wk',k(i,0)-i?.(k), pl{t)\, (Bl) 
with iJ,(k)"i'"^ = (5„„^,(yk^, - K^^AE), 



{7^k,k'(^,o)7^(;:;^p^(^) 
7^k,k'(^,o)p^,(^)7^^;:^ +{-}^ 



,n.q~ 

X S{nAE + £k' — ^k) 



(B2) 



dtpm 



ep 



\i' .71, X.zL.qz 

X (5(±WA,k-k',g, + nAE + £k' - £k)e'F'*'^^.k-'<',,, 



A,k-k%g; 



7^k,k'(^,o)(i-p^,(^))7^^;:)kP^(^) 



7^k,k'(^,o)p£,(^)7^^7Ul-Pk 



m} 



,(B3) 



5*pf(i)L, = - E ^[E^k-k',,J/(*g.)|^ 

k'.k", n,n' qz 

X S{nAE + Ek' - Ek + Ek" - ek"-k+k') 

X {Wk,k'(i,0)(i-p^,(t))wg2p^W 

X Tr [W^":k"'lk+k'P^"-k+k' (i) Wk"-k+k',k" (i, 0) 

X {i-pl,it))] - Wk.k'(i,0)p^,(i)w(?2 

X (i - P^W)TrK— - pf„_k+k'(t)) 



xWk"-k+k',k"(t,0)p^„(t)] } + {...}' 



(B4) 



Now, the kinetic spin Bloch equations can be treated 
via the numerical scheme in Ref. 42. The only differ- 
ence is the summations over sideband indexes which in- 
crease the complexity of the calculation. Typically, the 
sideband index runs through [—24,24] ([—3,3]) for the 
electron-impurity and electron-phonon scattering (the 
electron-electron Coulomb scattering) to converge the 
results when the THz field is E — 1.5 kV/cm with 

= 0.65 THz. 
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As mentioned in Sec. Ill, the initial distribution of the 
electron system is chosen to be a spin polarized hot- 
electron distribution under the THz field, which is ob- 
tained by sufficient long time (typically ~ 10 ps) evolu- 
tion from a spin-polarized Fermi distribution at the lat- 
tice temperature with the SOC being turned oS,*^ In Fig. 

we plot the time evolution of the distribution on the 
two Floquet states with k = (l.lfc^,0) {kp is the Fermi 

wave vector): p^^ •* (solid curve) and p^*^^^^ (dotted 
curve) for the case with 4 % initial spin polarization along 
the z axis. 



0.6 




t (ps) 

FIG. 14; (Color online) Time evolution of the distribution of 
two Floquet states with k — (l.lfc^,0): p^' ■* (solid curve) 
and Pk'^^' (dotted curve), (a): initial distribution prepa- 
ration without the SOC; (b): after the preparation with the 
SOC. E = 1.5 kV/cm, T = 50 K and N, = 0.05iVe. 

It is seen from Fig. [TiT a) that after only about 3 ps, 
the distributions show regular oscillations. This indicates 
that the system reaches its steady state. Moreover, the 
periods of the oscillations are close to the period of the 
THz field Tq. As we have pointed out in Sec. II, the eigen- 
modes of the steady-state distributions have the general 
form pg = e^'^k* g^."e™"* (a = 1,2,3,4) according 
to the Floquet theorenij^i The two diagonal elements of 
the distribution function should be governed by one of 
these modes while other modes are damping modes which 
do not appear in the steady state. When the THz field is 
not too large, the eigen- values of the relevant eigen-modes 
pjj are close to zero. Thus the distribution functions still 
have good periodic behavior and the period is close to 
To. 

From Fig. [HTbl one further finds that when the SOC 
is included, the steady-state distributions still have good 
periodic behavior. The system approaches steady state 



0.6 ■ 

0.4 ■ \ 

^ \ 

0.2 ■ X. 

^ TT'Tv-r^ 

20 40 60 80 

8 (meV) 

FIG. 15: (Color online) The hot-electron distribution in en- 
ergy space at t = 23 ps (solid curve) and the fitting curve 
(dashed curve). E = 1.5 kV/cm, T = 50 K, and Ni = 0.05iVe. 




0.2 0.4 0.6 0.8 1 



t-to ( To ) 

FIG. 16: (Color online) The hot-electron temperature Te 
(solid curves) and the chemical potential fj, (dashed curves) 
as function oit. E^ 1.5 kV/cm, T = 50 K and Ni = OmN^. 
Note that the scale for the dashed curve is on the right hand 
side of the frame. 

within 3 ps and the period is again close to Tq. The distri- 
bution difference on the two Floquet states in Fig. [HTal 
is due to the spin polarization whereas in Fig. IMf b) 
is caused by the spectral difference of the two Floquet 
states. 



APPENDIX C: HOT-ELECTRON EFFECT AND 
HOT-ELECTRON TEMPERATURE 

As has been shown in Appendix B that the steady-state 
distribution p^ (t) is a time dependent function which still 
exhibits good periodicity in our parameter regime. We 
can extract the distribution in energy space at any time 
t via Fourier transformation: 

F{cu,t)D{uj,t) = ^p^('"')(i)I)k„(^,i)- (CI) 

k?7 
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Here D{uj,t) = J^kri-^^vi^^^) the generalized density 
of states (2x2 matrix) where t is the center-of-mass 
time> 



.14.30 



/°° rlr T T 



(C2) 



It has been found that Dkri{uj,t) and D{Lo,t) are pe- 
riodic functions of t with the same period as that of 
the THz field Tq = 27r/fi. Therefore, the distribu- 
tion F{uj,t) is also a periodic function of t with pe- 
riod Tq. According to the symmetry analysis, D{uj,t) — 
D-ff(u;,t)l — Im{£>||(a;, i)}(Tj^.^ As the matrices i and 
ay form a group, the distribution F{uj,t) should also 
be decomposed into two parts: F{uj,t) = /(aj,i)l -I- 
s{uj,t)ay. Equation (|Cip then turns into [denoting 

f{LJ,t)D^^{u;,t) - s{u;,t)lm{Dn{io,t)} 

= C^H^,t), (C3) 

- s{u;, t)D^^{u, t) + f{u, t)lYn{D^^{uj, t)} 

= Im{CTH^,t)}. (C4) 

The solutions of the above equations are given by 

C"i?„-Im{CT^}Im{i^u} 



/ 



- [Im{i?u}] 
-Im{CTi}i?TT + C"lm{i?Ti} 



(C5) 
(C6) 



One notices that Eq. (CI) is a natural generaliza- 
tion of the distribution in energy space from thermal 
equilibrium to the nonequilibrium case. It is straight- 
forward to see that in the zero THz field limit the 
distribution F(uj,t) recovers the Fermi distribution as 

Dk,^{uj,t) = (5(w-ek^)|^kr,)('^k,7l and pl^'^'^\t) ^ fpiskri) 
in the zero- field limit {fp{x) is the Fermi distribution 
function). Therefore from Eq. (CI), F(uj,t) = /i^(w)i. 



A typical f{uj,t) is plotted in Fig. [121 We use the 
hot-electron temperature Tg to measure the hot-electron 
effect. The hot-electron temperature is determined by fit- 
ting the tail of f{Lu,t) with Fermi distribution function. 
We plot the fitted hot-electron temperature and the 
chemical potential ^ in Fig. [16] {to in the figure denotes 
the starting time which is 21.5 ps). It is seen in Fig. [TBI 
that Te and fj. are also periodic functions of t with pe- 
riodicity of To/2. This is because these quantities only 
depend on the strength of the THz field. The hot-electron 
temperature used in Sec. HI is the largest temperature, 
which is sufficient in measuring the hot-electron effect. 
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